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Abstract 

We want to understand he concentration of damage in microfractured elastic 
media. Due to the different scallings of the volume and area (or area and length in 
two dimensions) the traditional method of homogenization using periodic arrays of 
cells seems to fail when applied to the Mumford-Shah functional and to periodically 
fractured domains. 

In the present paper we are departing from traditional homogenization. The 
main result implies the use of Mumford-Shah energies and leads to an explanation 
of the observed concentration of damage in microfractured elastic bodies. 

1 Introduction 

A new direction of research in brittle fracture mechanics begins with the article of 
Mumford & Shah [15] regarding the problem of image segmentation. This problem, 
which consists in finding the set of edges of a picture and constructing a smoothed 
version of that picture, it turns to be intimately related to the problem of brittle crack 
evolution. In the before mentioned article Mumford and Shah propose the following 
variational approach to the problem of image segmentation: let g : C K 2 —>■ [0, 1] be 
the original picture, given as a distribution of grey levels (1 is white and is black), let 
u : — > R be the smoothed picture and K be the set of edges. K represents the set 
where u has jumps, i.e. u € C l (£l \ K, R). The pair formed by the smoothed picture u 
and the set of edges K minimizes then the functional: 

I(u, K) = [ a |Vu| 2 dx + f (3 | u - g | 2 dx + ^(K) . 

The parameter a controls the smoothness of the new picture u, 13 controls the L 2 
distance between the smoothed picture and the original one and 7 controls the total 
length of the edges given by this variational method. The authors remark that for (3 = 
the functional I might be useful for an energetic treatment of fracture mechanics. 
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An energetic approach to fracture mechanics is naturally suited to explain brittle 
crack appearance under imposed boundary displacements. The idea is presented in the 
followings. 

The state of a brittle body is described by a pair displacement-crack, (u, K) is such 
a pair if K is a crack — seen as a surface — which appears in the body and u is a 
displacement of the broken body under the imposed boundary displacement, i.e. u is 
continuous in the exterior of the surface K and u equals the imposed displacement uo 
on the exterior boundary of the body. 

Let us suppose that the total energy of the body is a Mumford-Shah functional of 
the form: 



The first term of the functional E represents the elastic energy of the body with the 
displacement u. The second term represents the energy consumed to produce the crack 
K in the body, with the boundary displacement uo as parameter. Then the crack that 
appears is supposed to be the second term of the pair (u, K) which minimizes the total 
energy E. 

Models for brittle damage, based on functionals of the Mumford-Shah type have 
have been proposed by Francfort-Marigo [TT], Buliga [6], among others. Such models 
have been studied intensively from the mathematical point of view, espacially by the 
Italian school of geometric measure theory, to name a few: De Giorgi, Ambrosio, Dal 
Maso, Buttazzo. 

The first homogenization result, concerning the Mumford-Shah functional, seems 
to be Braides, Defranceschi, Vitali [5]. In this paper it is done the homogenization of 
a Mumford-Shah functional of the form: 



The paper Focardi, Gelli |14j (and the references therein) are part of another line of 
research which might be relevant for this paper: homogenization of perforated domains. 

In the present paper we are departing from traditional homogenization. The line of 
research concerning perforated domains is close to our problem, but for various reasons 
the results from perforated domains don't apply here. 

We want to understand he concentration of damage in microfractured elastic media. 
Due to the different scallings of the volume and area (or area and length in two dimen- 
sions) the traditional method of homogenization using periodic arrays of cells seems 
to fail when applied to the Mumford-Shah functional and to periodically fractured 
domains. 

The main result, theorem 14.21 implies the use of Mumford-Shah energies and leads 
to an explanation of the observed concentration of damage in microfractured elastic 
bodies. 

Instead of performing a homogenization of the total energy of the microfractured 
body and then study the minimizers of the homogenized energy, we proceed along a 
different path. We study sequences of problems on fractured elastic bodies, indexed by 
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a scale parameter e. Each such problem has (at least approximative) solutions. We 
find estimates of the area of the damaged region in terms of the scale e. 



2 Notations 

Let Q be a bounded, open subset of R 2 , with locally Lipschitz boundary. We denote 
by Y = [0, l] 2 the unit closed square m K . 

For a given e > let Z £ C R 2 be the lattice of points in R 2 with coordinates of the 
form (em, en), for all m,n G Z. 

We denote by Z(e, 0) C Z £ the set of all z£Z £ such that 

z + eYc!] . 

To any z € Z(e, Q) we associate the cell 

D z = z + eY C Q . 

The set Z(e, Q) is finite for any e > 0. We denote the cardinal of this set by N(e) 
and we notice that as e goes to we have 

lim^ = l , 
e^o A(Q) 

where A(£l) denotes the area of 0. Thus for small e the number of cells N(e) is 
approximately equal to A(Q)/e 2 . 



3 The model 

We take Q to be the configuration set of a microfractured linear elastic body. We 
explain further what we mean by this. 

The elastic properties of the body are described by an elastic potential 

w : M 2 * 2 (R) -> R . 

We suppose that the function w is quadratic and strictly positive definite. 

For a given displacement u : — > R 2 , the elastic energy of the body is given by 

w(e(u)) dx , 

where e(u) is the deformation of the displacement u, that is the symmetric part of the 
gradient of u: for any x£l! 

e(u)(x) = ^(Vu(x) + (Vu) T (x)) . 

For a fixed e > we suppose that the body contains a distribution of micro-fractures 
at the scale e, seen as a union of (Lipschitz) curves 

F £ = (j (z + eF z ) , 
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where for each z € Z(e, £1) the (Lipschitz) curve F z lies inside the unit cell Y: 

F z C(0,lf . 

We explain further what we mean by an imposed boundary displacement uo, and 
what we mean by u = Uq on the boundary of f2. 

We consider, for simplicity, that uo : d£l — > M n is a continuous and therefore 
bounded function. Then, for any u € SBD(O), u = uq if the approximate limit of u 
equals uo in any point of 9f2 where the first exists, i.e.: for all x € d£l, if there exists 
v(x) such that 

,. $BJx)M I -v(x) I dy 
km — : — : = 

p^o+ I Bp{x) n ft 

then v(x) = uo(x). 

Definition 3.1 The class of admissible displacements with respect to the distribution 
of cracks F e and with respect to the imposed displacement uq is defined as the collection 
of all u G SBD(Q) such that 

(a) u = uo on d£l, 

(b) F £ CS U . 

This class of admissible displacements is denoted by Adm(F e , Uq). 

This definition deserves an explanation. An admissible displacement u is a function 
which has to be equal to the imposed displacement on the boundary of (condition 
(a)). Any such function u is a special function with bounded deformation, that is a 
reasonably smooth function on the set £l\S u and the function u is allowed to have jumps 
along the set S u . For the technical details see the Appendix. We have to think about S u 
as being a collection of curves, with finite length. Physically the set S u represents the 
collection of all cracks in the body under the displacement u. The condition (b) tells 
us that the collection of all cracks associated to an admissible displacement u contains 
F £ , at least. 

Definition 3.2 With the notations from definition \3.1\ the total energy of an admis- 
sible displacement u € Adm(F e , Uq) is given by 

E e (u)= [ w(e(u)) dx + GH\S U \F £ ) . 
Jn 

The energy of an admissible displacement is of Mumford-Shah type. It contains 
two terms. 

The first term measures the elastic energy of the body under the displacement u. 
Notice that in the expression of the elastic energy we have integrated over the whole 
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domain f2. This is simply because the collection of cracks associated to u (that is the 
set S u ) has Lebesque measure 0, therefore we have 



In physical terms, the right hand side expression would make more sense than the left 
hand side, but from the mathematical point of view they are the same. This is not 
meaning that the elastic energy neglects the fractures. Indeed, further we shall infimize 
the energy E £ over the whole set of admissible displacements. According to condition 
(b) of definition 13 . 1 (. this set is defined with respect to the collection of cracks F £ , 
therefore the infimum of the energy E £ depends on the set of cracks F £ . 

The second term of the Mumford-Shah energy measures the surface energy caused 
by the apparition of new cracks. The collection of new cracks is the set S u \F e . The 
constant G has the dimension of energy per unit area, and it is physically related to 
the Griffith constant. 

In [3] has been proven that functionals like E £ are L 1 inferior semi-continuous 
and coercive, hence on closed subspaces V of SBD(S7) the functional E £ has a mini- 
mizer. Such a closed subspace of SBD(f2) is the space of all admissible displacements 
Adm(F £ , uq) . Therefore we have: 

Theorem 3.3 On the space Adm(F £ , Uq) we consider the topology given by the con- 
vergence: Uh — > u if 



Then there exists a minimizer of the functional E £ over the set Adm(F £ , uq). 

In the following section we shall use approximate minimizer s. 

Definition 3.4 For a given 5 > 0, a function u € Adm(F £ , uq) is a S- approximate 
minimizer if 



For fixed S > 0, we model an approximate displacement of a microfractured body 
as a sequence of displacements u e , with e converging to 0, such that for each e > the 
displacement u e € Adm(F £ , uq) is a <5-approximate minimizer of the Mumford-Shah 
energy E £ , over the set Adm(F £ , uq). 

Notice that in the model, at this stage, there is no relation between the crack sets 
F £ ,F £ i, for two different scales e,e'. 

4 An estimate related to damage concentration 

For fixed e, 5 > 0, given F £ and imposed boundary displacement uo, let u € Adm(F £ , uq) 
be a 5-approximate minimizer of the Mumford-Shah energy E £ . 





E e (u) < S + M{E £ (v) : VG Adm(F s ,u )} 
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In this section we want to estimate the number of e-cells z + eY , z 6 Z(e, $7), where 
the initial cracks z + ei 7 ^ propagated. 
Let I > be a given length. 

Definition 4.1 For any ceZZ Z? z = z + eY, z £ Z(e, f&), and any 5 -approximate mini- 
mizer u we define the emergent crack in the cell D z by 

S u (z) = (z + eY)n(S u \(z + eF z )) . 

A cell D z is called active if the length of the emergent crack is greater than el, that is: 

n\S u {z))>el . 

We denote by M(e,l) the number of active cells. (In this notation we don't mention 
the dependence of M(e, /) on the 5-approximate minimizer u.) 

Theorem 4.2 Suppose that for fixed 5 > 0, the crack sets F e are chosen so that there 
exists an approximate displacement of a microfractured body u E , with e converging to 
0, with the property that the sequence 

inf{£ £ (v) ; vG^dm(F £ ,u )} 

is bounded. 

Then the number of active cells M(e, I) is of order 1/e and the area of the damaged 
region of the body 

Damaged(e, Q) = [^J D z 
D z active 

is of order e. 

Proof. Let M > such that for all e > we have 

inf {E £ (v) : v € Adm(F £ ,u )} < M . 
According to definition 13.41 for any e > we have 

E £ (u £ ) = [ w(e(u £ )) dx + GH 1 (S Ue \ F £ ) < 
Jn 

< o + inf {E e (y) : v G Adm{F £ ,u )} < 5 + M . 
From definition 14.11 we get the following estimate: 

H 1 (S Ue \F £ )= H 1 (S u (z))) > M(s,l) I e . 

We have therefore 

G M(e, I) I e< GH 1 {S Ue \ F £ ) < E £ (u £ ) <M + 5 . 
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All in all we have obtained the estimate: 

n 1M + S 

The area of the damaged region of the body is 

Area(Damaged(e,£l)) = ^ Area(D z ) = e 2 M(e, I) < ~ ■ 

D z active 

The proof is done. □ 

5 Conclusions 

The theorem implies that the area of the damaged region is much smaller than the total 
area of the body, as e goes to zero. In this model the use of Mumford-Shah energies 
leads to an explanation of the observed concentration of damage in microfractured 
elastic bodies. 

Notice that we need more precise estimates in order to prove that the damaged 
region (at the scale e) converges, as e goes to zero, to a curve with finite length. All we 
know at this moment is that the area of the damaged region goes to zero as the scale 
parameter e. 

In experiments it has been observed that the damaged region is approximately 
straight. It is possible that Mumford-Shah energies might explain this, since geometries 
of the active crack set, that is S U \F £ , close to a straight line would be preferred by the 
energy E £ . See [7] for examples that in some situations the leading term of a Mumford- 
Shah energy is the one accounting for the length of the crack, and not the elastic energy 
part. 

Finally, in theorem 14.21 we obtained an estimate of the number of cells where cracks 
of length at least e I appear. It would be interested to study the interplay between e 
and I in this estimate. 



6 Appendix. Functions with bounded variation or defor- 
mation 

This section is dedicated to a brief voyage trough the spaces SBV and SBD. 

The space SBV(f2, R n ) of special functions with bounded variation was introduced 
by De Giorgi and Ambrosio in the study of a class of free discontinuity problems ( [9] , pQ , 
[2]). For any function u € R n ) let us denote by Du the distributional derivative 

of u seen as a vector measure. The variation of Du is a scalar measure defined like 
this: for any Borel measurable subset B of f2 the variation of -Du over B is 

( oo 

I Du I (B) = sup < I Du{Ai) | : A.CB, A4 H Aj = Vi + j 

[i=i 
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A function u has bounded variation if the total variation of Du is finite. We send the 
reader to the book of Evans & Gariepy [13] for basic properties of such functions. 
The space SBV(f2, R n ) is defined as follows: 

SBV(ft, R n ) = {u G L^O, R n ) : \ Du | (Si) < +00 , | D s u | (SI \ S u ) = 0} . 

The Lebesgue set of u is the set of points where u has approximate limit. The com- 
plementary set is a C n negligible set denoted by Su- If u is a special function with 
bounded variation then Su is also a (i.e. countably) rectifiable. 

From the Calderon & Zygmund [8] decomposition theorem we obtain the following 
expression of Du, the distributional derivative of u G SBV(r2,i? n ), seen as a measure: 

Du = Vu(x) dx + [u] ® n dH^ 1 . 

We shall use further the notation A if the measure \x is absolutely continuous 
with respect to the measure A. 

Let us define the following Sobolev space associated to the crack set K (see [3]): 

W 1 / = ju G SBV(Sl,R n ) : J I Vu | 2 dx + J [u] 2 dH^ 1 < +00 , | D s u |< H'^ 1 

It has been proved in [TO] the following equality: 

w 1 ' 2 (n\K,R n ) nL°°(n,R n ) = w^ 2 (n,R n )nL°°{n,R n ) . (6.0.1) 

A similar description can be made for the space of special functions with bounded 
deformation SBD(O) can be found in [4j. For any function u G L 1 (J7, R n ) we denote by 
Eu the symmetric part of the distributional derivative of u, seen as a vector measure. 
We denote also by J u the subset of f2 where u has different approximate limits with 
respect to a point-dependent direction. The difference between Su and Ju is subtle. 
Let us quote only the fact that for a function u G SBV(f2, R n ) the difference of these 
sets is H n_1 -negligible. 

The definition of SBD(O) is the following: 

SBD(ft, R n ) = {u G L\n, R n ) : \ Eu \ (SI) < +00 , \ E s u \ {SI \ J u ) = 0} . 

If u is a special function with bounded deformation then J u is countably rectifiable. 
We have a decomposition theorem for SBD functions, similar to Calderon & Zygmund 
result applied for SBV functions. The decomposition theorem is due to Belletini, 
Coscia & Dal Maso [1] and asserts that 

Eu = e(u)(x) dx + [uJOndHp 1 . 

Here means the symmetric part of tensor product and e(u) is the approximate sym- 
metric gradient, hence the approximate limit of the symmetric part of the gradient of 
u. 

We sum up the main facts about functions with bounded variation or deformation, 
in the following three theorems. 
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Theorem 6.1 Let u G L 1 (p,,W m ). Then 

- (De Giorgi) Ifu G BV(ft,IR m ) then S u is countably rectifiable, H n_1 (S u \Ju) = 
and in H 71 ^ 1 -almost every point x G S u exists the approximate limits of u in 
the directions v{x) and —u(x) where v{x) is the normal to Su in x. 

- (Kohn, Ambrosio, Coscia, Dal Maso) Let m = n and u G BD(f2). Let O u be the 
Kohn set : 

@ u = Len: limsup |£»l W*)) > ) 

T/ien U is countably rectifiable , J u Q ©u 7i n_1 (0 u \ Ju) = . 
Theorem 6.2 Let u G L 1 (J7,R m ). Then 

- (Calderon, Zygmund) If u G BV^IR" 1 ) £/ien u is approximately differ entiable 
C n -a.e. in fi. The approximate differential map x \— > Vu(i) zs integrable. Du 
splits into three mutually singular measures on f2 

Du = Vu dx + [u] <g) z/ft! 1-1 + Cu 

where [u] zs the jump ofu in respect with the normal direction on S u Cu is the 
Cantor part of Du defined by Cu(A) = D s u(A \ S u ) where D s u is the singular 
part of Du in respect to C n . 

- (Belletini, Coscia, Dal Maso) Let m = n and u G BD(S7). Then u has symmetric 
approximate differential e(u) C n -a.e. in Q and Eu splits into three mutually 
singular measures on U 

Eu = e(u) dx + [u] vH^ 1 + E c u 

Ju 

Moreover u is approximately differ entiable C n -a.e. in f2. 
Theorem 6.3 The following are true: 

- W 1,1 ^, W 71 ) C BV(f2,R m ). The inclusion is continuous in respect with the 
Banach space topologies. If 

u G SBV(n,M m ) 

then 

uGW 1 ' 1 (0\S u ,R m ) 

Moreover if u G W 1 ' 1 ^ \ K,R m ) n L°°(n,R m ) , where K is a closed , countably 
rectifiable set withH n ~ l {K) < +oo, then u G SBV(f),IR m ) and H n ~ l {K \ S u ) = 
0. 

- Let LE l (£l) be the Banach space of L l (£l,W l ) functions with L 1 symmetric dif- 
ferential. If u G SBD(f2) then u G LE 1 (Q, \ J u )- Let K be a closed , count- 
ably rectifiable set with U n - l {K) < +oo. If u e LE 1 (ft \K)n L°°(n,R n ) then 
u G SBD(ft) and H n - l {K \ J u ) = 0. 
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